Non-Gaussianity from two right-handed sneutrino curvaton decays 
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In this paper, we consider the effect of two right-handed sneutrino curvaton decays and investigate 
the parameter space. We compare the difference of the result between single- and two-curvaton cases. 
We find one Yukawa coupling of the right-handed snetrinos can be as large as A ~ 0.1 while the 
other one is much smaller. When the curvatons decay, we assume both of them subdominate the 
energy density of the universe. We find that, unlike a single curvaton case, here a small or negative, 
as well as a large fNL can be generated. 



CN ' PACS numbers: 

o 

(N 
513 



I. INTRODUCTION 



■ Although the idea of inflation [ll-Q (for a review, see Q ) started from solving the problems (like horizon and flatness 
OA ' problem) of the old hot big bang model, the key to distinguish different inflation models lies on primordial curvature 
perturbation C. which provides the seeds for structure formation. Inflation could more or less generate some primordial 
' curvature perturbation since we are living in a quantum world. The primordial curvature perturbation generated 
rS I from simple (single-field slow-roll) inflation models is adiabatic, almost scale invariant, and Gaussian. However, if 
JL ' inflaton is responsible for generating primordial curvature perturbation, the constraint of CMB normalization (i.e. 

1/2 

. P^' ~ 5 X 10-^) is so strong that we usually need to fine-tune the parameter(s) when building an inflation model. 
' Furthermore, future experiments (like PLANCK satellite) may detect large non-Gaussianity and hence rule out simple 
scenarios of inflation. This will reveal for us a nontrivial way of generating primordial curvature perturbation and 
' one of the promising ideas to generate large non-Gaussianity is through a curvaton Non-Gaussianity generated 

^ , from curvaton scenario can be described by the nonlinear parameter /nl, which takes the form 
OO 

c=c, + lfNLe, + ---, (1) 

jy-^ ' where Q denotes the Gaussian part of C,. Currently the upper bound of f^L is roughly given by (2(t) [8| 
^ ■ \fNL\ < 100. (2) 

In the near future, the PLANCK satellite will reduce the bound to IJnlI < 5 if non-Gaussainity is not detected. 

^ Curvaton is supposed to be light^ and subdominant during inflation. Because it is light, it can produce sizable 

' quantum fluctuations which when stretched outside the horizon during inflation would become classical perturbations. 
$^ . Because it is subdominat, the perturbations should be regarded as isocurvature perturbation after inflation when the 
curvaton field starts to oscillate. The curvaton is supposed to decay after inflaton decay and at the same time 
transform its isocurvature perturbation to curvature perturbation. If the curvature perturbation of the universe is 
from the curvaton, it could liberate the constraint on inflation and lower the scale of inflation It would be 
interesting if we can identify a field from particle physics to be the curvaton and one of the candidate is a right- 
handed (RH) sneutrino in the framework of supersymmetry (SUSY). The idea of a right-handed sneutrino curvaton 
has already be mentioned in the original paper of curvaton [5[ and been considered in In this case, however, 

it is possible that we have more than one generation of the right-handed sneutrino and more than one curvaton. The 
calculation for two curvaton decay has been considered in [18|, however, the parameters used in the paper is not of 
direct use for particle physicists to find possible candidate for the curvatons. In this paper we consider the decay 
width of the right-handed sneutrino curvaton in terms of the Yukawa couplings and masses. It turns out that there 
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are six parameters in addition to the Hubble parameter. We then numcrieally scanning the parameter spaee in order 
to find solutions. 

This paper is organized as foUowsing. In section we review the result of [3| in order to compare single- and 
two- curvaton decays in the succeeding sections. In section Hill we present the formalism we use for calculating two 
curvaton decay. In section ITVl we specify the parameter space and present our numerical results. Section |V] is our 
conclusion. For completeness, there is also an Appendix section in the end of the paper where we summarize the 
equations adopted from [l8| . 



II. SINGLE RIGHT-HANDED SNEUTRINO CURVATON 



It is well known that in standard model of particle physics, there is no good candidate to play the role of a curvaton. 
However, if we go beyond the standard model by imposing SUSY, there are lots of scalar fields. In this section we 
review the idea of a right-handed sneutrino which plays the role of the curvaton. We focus on the case that when 
the curvaton decays, it subdominates the energy density of the universe [l6j . because large non-Gaussianity could be 
produced in this region of the parameters. 

The superpotential of the RH neutrino is given by 



m<i>2 



(3) 



where 3> is the RH neutrino superfield, and L are the MSSM Higgs and lepton doublet superfields, and m is the 
RH neutrino mass. The canonical type-I seesaw mechanism gives the mass relation between light and heavy Majorana 
neutrino masses as 



(4) 



lO^GeV denoting the vacuum expectation value of Hu- 



with Vu 

neutrino oscillation data, Amf2 



7.59t° l^i X 10- 



^eV^ and |Am| 
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The mass squared differences revealed from the 
= 2.43 ± 0.13 X lO-^eV^ [ll, indicates that if we 



consider three generations of RH sneutrinos, the lightest left-handed neutrino mass can be very small while the mass 
of the other two generations are fixed due to the oscillation data. This means in principle we can have a Yukawa 
coupling arbitrarily small (7713 ^ 0) but the other two Yukawa couplings should result in the neutrino masses around 
7711^2 ^ O.leV (with a mass difference Ami2 ^ 10~^ eV, namely, "inverted ordering"). On the other hand, we can 
also have "normal ording" which implies mi ~ 0, m2 ~ 0.01 eV, and TO3 ^0.1 eV. From Eq. ([4]) we can see that for 
the case - 0.1 eV, if we have m - lO'^Mp (to ~ lO'^Mp), we need - 10"^ (A^. - lO'^). 
The potential of the RH sneutrino a can be expressed as 



Via) 



1 



(5) 



For simplicity, here we do not consider the Hubble-induced mass term. It is possible that the mechanism which 
suppresses the Hubble-induced mass term for the inflaton also suppresses that for the curvaton. For example, this is 
the case for D-term hybrid inflation [l^. The decay rate for the RH sneutrino is 



The spectrum is given by 



1/2 



1 = TO. 

4tt 



5 X 10" 



(6) 



(7) 



If we assume that at the time to of curvaton oscillation with energy density Pa{to) = "^^o'»/2, the universe is 
dominated by radiation (the decay products of inflaton) with energy density PR{to) = im^Mp. At the time of 
curvaton decay to, the energy density of the universe is given by PRitn) = 3r^M|, = PR{to){a{to) / a{tD))'^ ■ Therefore 
a{tD) /a{to) = (to/F)^/^ and ^a,D is given by 



Ptot 



6 \Mp J \r 



6 \Mp J K 



(8) 
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Throughout this paper, we always use a subscript to denote horizon exit. When Vla-.D is small, the nonlinear 
parameter is given by 

Note that here /nl can never be negative because we assume the eurvaton is subdominant when it decays. As we will 
see in the following section, this is not true in the two-curvaton case. By using Eqs. (O, (HJ, and (|9]), we can obtain 

X^M^ = 3.9 X lO^cr,^, (10) 

and 

fNL^ (8.25x10^)^. (11) 

It is interesting to note that here the spectrum does not give constriant on the mass of the right-handed sneutrino. 
However, as we can see in the next section, the results does depend on the masses when we consider the two curvatons 
case. 



III. TWO RIGHT-HANDED SNEUTRINO CURVATONS 



In our setup, we consider three generations of right-handed sneutrinos for type-I seesaw model and for simplicity 
we assume the heavist right-handed sneutrino (with a mass we denoted as nic) does not play any role in cosmology, 
namely rric > H^? Here we consider the case that the lighter two right-handed sneutrinos are curvatons (we call 
them eurvaton a and eurvaton h) and investigate the effects on the primordial curvature perturbation. 

The decay rates of the two curvatons are 

Ta^-r^a and rfa = -Amb (12) 
47r 47r 

where Xa and A;, are the Yukawa couplings and ma and ruh are their masses. Quite generally, we consider 

H^> ma> mb>Ta> Fb. (13) 

We assume that when both of the eurvaton decays, the energy density of the universe is dominated by radiation so 
that we can compare this with single eurvaton case. The similar calculation of the energy density ratio SI as in Eq. ([8]) 
can be obtained for radiation 7, eurvaton a, and eurvaton h. At eurvaton a decay (which is denoted by subscript "1" 
throughout this paper). 



^al — 



1/ \2/ \l/2 1/ \2 rr- 



6 \MpJ yVaJ 6 \MpJ Xa ' 

1 / J, \ 2 / \ 1/2 1 / I, \ 2 / < \ 1/2 



where the subscript 70 denotes pre-existing radiation just before eurvaton a decay. At eurvaton b decay (which is 
denoted by subscript "2" throughout this paper), 

6 Mp Xa 



^ If the gauge non-singlct scalar fields like a Higgs arc light during inflation, the fields can establish an expectation value through quantum 
fluctuation. Through a Yukawa coupling, this gives an effective mass to the right-handed sneutrino field as can be seen in Eq. 151 . 
The effect is most significant for the heavist right-handed snetrino due to the largest Yukawa coupling. Even if m,c < H, it is possible 
that the heavist right-handed sneutrino obtains an effective mass larger than the Ifubblc scale during inflation while other right-handed 
sneutrinos have masses lighter than the Hubble scale. 
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where the subscript 71 denotes radiation just before curvaton b decay. 
To linear order, the spectrum is given by [l8| 

^(1) = ^'^Cu, + 

where the parameters A and B can be found in the Appendix and 



1/2 



1 H, 



and P, 



1/2 



1 



(16) 



(17) 



Here we use subscript "(1)" to denote "first order" and "(2)" will be used to denote "second order". If we define 

/3 = a,/6,, we obtain 



We can write the spectrum explicitly by inserting A and B to obtain 

P^[i-j = P^^^-^{H^^.,a^,b^,ma,nib,Xa,\b) 



(18) 



3^ 



\ 



12M2 0Ff^ + 



Stt, 



rria ^ A„ A A„ 



'nib. _L JLA ^ / '»b 



4A4 + 0F(^ + ,/^^) 



Aa 



ma Aa 



2 r 



m/|, + v^(4^ + 3^) 



(19) 



This is subjected to CMB normalization P^l^i^ ~ 5 x 10 "''at horizon exit. 
To second order, the curvature perturbation is given by [isj 

C = C2 = C2(l) + ^C2(2) = [^Ca(l) + -SCh(l)] + ^ ^Ca(l) + -^Ch(l) + £'Ca(l)C6(l) 



(20) 



where ^2 is the total curvature perturbation after the second curvaton decay. As in Eq. (|17p . we use subscript "(1)" to 
denote "first order" and "(2)" to denote "second order". The first order part (2(1) is Gaussian and the second order 
part (^2(2) is non-Gaussian. The parameters C, I?, E can be found in the Appendix and the nonlinear parameter is 



5 + ^I3^EAB + li^DB'^ 
" 6 (A2 + /?2B2)2 



(21) 



IV. NUMERICAL RESULTS 

There are six parameters: a^:,bt,,ma,mb, Xa, in addition to the Hubble parameter ff* and one constraint (CMB 
normalization). We tackle the problem numerically by scanning the parameter space. First of all, we assume the 
inflaton contribution to the curvature perturbation is small, this implies 



P, 



1 



2V27r JenMp 



< 5 X IQ- 



(22) 



where en = —H/H^. For a typical value of en ~ 0.01 at horizon exit, we have the bound for Hubble parameter 

iJ* < V2tt X lO^^Afp. (23) 

We will numerically find solutions by making plots for the two cases iJ* = lO^^Mp and = lO^^Mp in this paper. 
The range of the curvaton field values at Hubble exit are chosen to be 

-ff* < a* < A/p and <b^ < Mp. (24) 

The lower bound is from the requirement that the classical field value is larger than its fluctuation and the upper 
bound is from the requirement that the curvaton would not drive a second stage of inflation. The masses of the 
curvatons are 



10 ^^Mp<ma<H^ and 10 '^'^Mp < mi, < rua 



15 J 



(25) 
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The lower bound is from the fact that the mass of the RH sneutrino should be larger than its soft mass which is 
assumed to be TeV scale and we also use the constraint from Eq. ([T^ . The decay rate is also constrained by Eq. 



therefore we choose 



And again we require 



therefore we choose 



47r 



< rub, 



10"'" < Aq < min{J— 47r,l}. 



)2 y2 

Tq = -^ma >Tb^ ir'^b-, 



(26) 



(27) 



(28) 



10" 




(29) 



We plot our results in Figs. [T][3] for = lO^^A/p and Figs. I311S] for = lO^^A/p. For comparison, we also plot the 
results from the single-curvaton case. Naively one may think it is possible that one of the curvaton plays no role at all 
and the results would be dominated by a single curvaton. However, those plots show that both of the curvatons could 
play some role in generating curvature perturbation. In all of the plots, each point represents a solution of all six 
parameters subjected to CMB normalization. We pick up a few points and list them on the TableU) In the plots Jml is 
an output after we imposing CMB normalization. We notice that even we assume the energy density of the curvatons 
is subdominant when they decays, we can still get small or negative (as well as large) /jvl- In addition, we also find 
that one of the Yukawa coupling can be as large as ~ 0.1 while the other coupling is relatively small. This is interesting 
because it is consistent with the neutrino oscillation data which we discussed below Eq. (|4|) by the assumption that 
the curvatons are RH sneutrinos. For example, in the second row of Table HI we would have the light neutrino masses 
correspond to the heavy neutrino masses ma and to;, as m^a ^ n^i ^0.1 cV and m,yb "m-a 10~^ eV respectively 
from Eq. Q. We can easily choose a mc > ma with a larger Yukawa coupling to make m^^c m2 ~ 0.11 cV and make 
an inveted ordering of neutrino masses compatible with neutrino oscillation. As another example, for the seventh 
row of Table HI m^a ~ '7i2 ^ 0.01 eV can be obtained and along similar arguing we can obtain a normal ordering 
of neutrino masses. We would like to emphasis here that our goal is NOT showing that all of our parameter spaces 
are compatible with neutrino oscillation data, because our results can generically apply to other two-curvaton models 
with a similar decay rates. However, it is interesting enough if some of them do compatible with our assumption that 
two RH sneutrinos can play the role of curvatons. 
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TABLE I: List of a few points of the numerical solutions. 
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FIG. 1: Xa versus Ai,. The points outside the range of single-curvaton bound represent the effect of considering two curvatons. 
As we can see here, one of the Yukawa couphngs can be quite large. 
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FIG. 2: a, versus fe*. 



V. CONCLUSION AND DISCUSSION 



We have explored the parameter space of generating primordial curvature perturbation via two right-handed sneu- 
trino curvaton decays in the framework that there are three generations of RH sneutrinos. We compared the results 
with the single-curvaton case and found that an additional RH sneutrino curvaton could have some effects and cannot 
be neglected. Notably we may still get small or even negative Jn l in the case that the energy density of the curvatons 
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FIG. 4: Xa versus At. We found fewer points because it is more difficult to find solutions by our numerical computation for 
smaller Jf*. 



sub-dominates when they decay. 

As can be seen from Fig. [5] and Fig. [5] that the field value for the solutions may be not very small compared to the 
Planck mass (although it has to be smaller than the Planck mass in order to avoid driving a second stage of inflation) . 
Therefore we have to suppress the nonrenormalizable terms in the superpotential. This can be achieved for example 
by judiciously assigning R-charge to the RH sneutrinos [lll[20|. 
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Appendix A: Notation 



The purpose of this appendix is to show how to obtain parameters A, B, C, D, and E from Eqs. and (|15l) . 
The following equations and notations are from [isj : 

The energy density parameters (fli = pi/ptot)- 

flai the density of the first curvaton a just before the first decay. 
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flbi the density of the second curvaton b at the first decay. 

^{,2 the density of the second curvaton b just before the second decay. 

51^0 1 the density of pre-existing radiation at the first (curvaton) decay. 

ri^ji the density of all radiation immediately after the first decay. 

r2^j2 the density of all radiation just before the second decay. 

r2^2 2 the density of all radiation immediately after the second decay. 

The full non-linear curvature perturbations: 

Ci— C2) the primordial perturbation after the second curvaton decay, but before the nucleosynthesis, 

Ci the total perturbation at the first decay, 

(2 the total perturbation at/after the second decay, 

the perturbation of the first curvaton a, 
Cf, the perturbation of the second curvaton b, 
C70 the pre-existing radiation perturbation, 

the radiation perturbation after the second decay. 

Subscript in parenthesis: 

the first order part of the primordial perturbation, 
C(2) the second order part of the primordial perturbation, 
Ci(i) the first order part of the total perturbation at the first decay, 
Ci(2) the second order part of the total perturbation at the first decay, 
C2(i) the first order part of the total perturbation at the second decay, 
^2(2) the second order part of the total perturbation at the second decay, 
Ca(i) the first order part of the first-curvaton perturbation, 
Ca(2) the second order part of the first-curvaton perturbation, 
Cb(i) the first order part of the second-curvaton perturbation, 
Cb(2) the second order part of the second-curvaton perturbation, 
C-yo{i) the first order part of the pre-existing radiation perturbation, 
Cfo(2) the second order part of the pre-existing radiation perturbation, 
Cfi{i) the first order part of the radiation perturbation after the first decay, 
Cyi{2) the second order part of the radiation perturbation after the first decay, 
the first order part of the radiation perturbation after the second decay, 
C-y2{2) the second order part of the radiation perturbation after the second decay. 

We have 

Cl(l) = /7olC7o(l) + ./alCa(l) + .fblCb{l) (Al) 

where 

^•^"^ ^ 4r!7„i + 3oll + 3r!w' ^^^^ 

f = f A3) 

^ 4r!7„i + 3r!al + 3^61 ' ^ ' 

4r27oi -|- 3fiai + 3f2bi ^ 

After the first curvaton decays, but before the second curvaton decays, the curvature perturbation in the radiation 
at first order is given by 

C7,(i)=i?iCi(i)-(i?i-l)Cfc(i), (A5) 

where 



3 + /al 



3(l-/bi) + /ai 



(A7) 
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Compare Eq. (|Il8l) with Eq. we conclude 
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The non-linearity parameter is: 



C ^ C-^F, (A24) 
D = D-^G (A25) 



6 {rl+P^rlY 



-J ' 2'a) I ^'a'O > H ' 2'"' f\OR\ 

JNL^-^ —^772 l^^oJ 
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